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ABSTRACT 


Asymptotic ray theory, as presented by Cerveny and 
Ravindra (10) or Hron and Kanasewich (16), provides a 
conceptually clear and simple theory to solve a variety 
of problems dealing with the propagation of head waves in 
common media. The use of the theory is limited, however, 
due to singular points in the spatial domain. With this 
limitation in mind, it is desirable to calculate a more 
exact solution which lends itself to easy numerical calcula- 
tion. The integrand in these expressions can be obtained 
by satisfying boundary conditions. By a suitable deforma- 
tion of the original integration contours, and use of the 
method of steepest descent, results are obtained. The 
formulae give a continuous amplitude at the critical point 
of the ray and tend to the results obtained from asymptotic 
ray theory for epicentral distances beyond the end of the 
zone where pure reflected and head waves interfere. 

Amplitude-distance curves for various interference 
waves will be presented. The separation of head waves and 
reflected waves in a synthetic seismogram near the criti- 


cal point will also be given. 
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CHAPTER 1 


INTRODUCTION TO HEAD WAVES 


In the time domain, asymptotic ray theory provides 
a solution to the wave equation and boundary conditions 
in the form of an asymptotic series in inverse powers of 
frequency. Details of this method can be found in Hron 
and Kanasewich (16) and Cerveny and Ravindra (10). The 
ray theory is useful in finding expressions for the 
reflected and head waves in many problems which prove 
difficult to solve by integral methods (such as curved and 
dipping sen E aces or inhomogeneous media). It will be 
shown in this thesis that for a situation in which both 
the integral solution and asymptotic ray solution are 
known, they turn out to be the same for high frequencies. 

The general equation of motion for waves in a homo- 


geneous, perfectly elastic and isotropic medium is 


aan 2 
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where W is the particle displacement vector and t is the 
time. The time-harmonic solution to this equation is 
assumed expressible in inverse powers of frequency accord- 


ing to asymptotic ray theory 
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W = exp(i w(t-t)) , (i ms Wy 1.2 
k=0 


where t and Wy are independent of t and w. 

This expansion is called a ray series, t a phase 
function, and Wy the amplitude coefficients of the series. 
Moving surfaces of constant phase, t = t(x,y,z) will be 
called wavefronts and the orthogonal projections of these 
surfaces will be termed rays. 

If the solution for equation 1.1 is known for harmonic 
waves, the results can easily be generalized for impulsive 
waves using Fourier Transforms. The impulsive solution, 


denoted by W can be expressed as 


where W is the steady state solution of 1.1 and S(w) is 
the spectrum of the impulse. From 1.2 and the definition 
of Ww we find that 


where 


F(t - t) = ( (i oy S(w) exp[Li w(t - t)] dw 
0 , 
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It can be seen from this definition of f, that 
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Cerveny and Ravindra (10) state that this ray series 
is invalid at points where the function t(x,y,z) is non- 
analytic. Therefore the ray series is invalid where the 
time-distance curve of the wave has discontinuous deriva- 
tives, at a triplication point, for example. The invalid 
point to be considered in this study is the critical point. 
At this point on the surface the travel-time curves of 
reflected and head waves are still continuous, although 
they are tangential to each other. This tangency implies 
that the critical point is a multiple valued point. Beyond 
the critical point the head wave retains its linear travel- 
time relation while the reflected wave still maintains a 
hyperbolic relation. It will be shown that the region of 
invalidity (called an interference zone or critical region) 
of the ray series depends on the source frequency. 

Cerveny and Ravindra (10) have used asymptotic ray 
theory to derive the following amplitude formulae for head 


and reflected waves: 
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where pegs velocity of propagation of the m-th ray 
segment 
on = angle of incidence of the m-th ray segment 


on the interface upon which it impinges 


h_ = thickness of the layer housing the m-th ray 


m segment 

k = total number of ray segments 

f = source frequency 

Yr = epicentral distance 

r* = critical distance for the ray. This is--the 


minimum distance at which the head wave can 
physically exist. It corresponds to the ray 
which is critically reflected. 
The product term in the numerator represents the product 
of all reflection, transmission, and head wave coefficients 
encountered by the ray during its propagation from source 
to receiver. It is obvious that the head wave amplitude 
Tends Lonimntrimnity asere> re. 

A theoretical analysis of the region near the critical 
point is of great practical seismological importance 
(Musgrave (21)). First, the local maxima of the amplitude 
of the reflected waves lies in the region of the critical 
points and therefore are at the most favorable position 
to detect these strong waves. Second, if the critical 


point can be determined by experimental data, then the mean 
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velocity of the overburden can be calculated by examining 
the travel time curves of the head waves. Finally, on the 
basis of kinematic data alone, a seismic interpreter can 
select a series of equally true earth models, but will not 
be able to choose the correct one. If the maximum of the 
amplitude distance curves can be found experimentally and 
the position of the critical point calculated for several 
earth models, then the correlation Macy Gs actual value 
and the calculated one can be used as proof of the validity 
of the earth model chosen. 

It is the aim of this thesis to conduct a theoreti- 
cal and numerical investigation of the reflected waves 
which emerge near critical points. Expressions for the 
dynamic properties of interfering reflected and head waves 
will be derived. 

Head waves have been given several different names 
in the literature. These include conical, lateral, 
refracted, and Mintrop waves. The study of head waves 
from a plane interface between two homogeneous media has 
been well done in the field of seismic theory, for example 
Muskat (22). Head waves were later studied using the 
Sommerfeld integral method. The most concise of these 
studies is Heelan (14). Studies of the interference 
reflected wave-head wave system in the critical region 
have been considerably rarer. This was first investigated 
by Brekhovskikh (4) for liquid media. Cerveny (6) studied 


this interference phenomena for solid media for a single- 
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layered model with an arbitrary velocity distribution at 
the refracting interface. The arbitrary velocity distri- 
bution allowed one to consider media in which either P or 
S type head waves could exist. Later, Cerveny (9) solved 
the problem for a multi-layered medium, but with only one 
critical point in existence. Thus, the problem of solving 
interference wave amplitudes for multi-layered media and 
arbitrary ray paths remains unsolved. (The problem of 
multiply-reflected head waves is also in need of investiga- 
tion.) This thesis aims to present a solution to these. 
This solution is for horizontally-stratified media with 
all layers homogeneous. The formal solution is in the form 
of a single integral with a simple integrand for the one 
interface problem and an extremely cumbersome one for 
nents] A ayietted.enidad a These integrals have been evaluated 
exactly by several authors for the one interface problem 
(see Cagniard (5), Dix (11), Bortfeld (2), (3), Muller 
(20)). An excellent bibliography and a simple picture of 
head waves are given in Cerveny and Ravindra (10). 

A head wave is basically formed by the concave para- 
bolic wave front impinging upon a plane interface whose 
refractive index is less than one. The ray must strike 


the interface at the critical angle 6” defined by 


. 7 * e .J 
sin 6 =n = refractive index. 


Two criteria necessary for head wave formation are: 
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(1) the radius of curvature of the incident wave 
front must be finite (plane waves cannot gene- 
rate head waves) 

(2) the refracted segment of the ray must propagate 
with the greatest velocity of the entire ray. 

Every head wave is closely associated with a reflected 

and transmitted wave. The head wave is a plane wave whose 
end points are the point where the transmitted wave touches 
the interface and the point of tangency of head wave and 
reflected wave. This is shown diagramatically in Fig. 1.1. 
This thesis concerns itself only with the head waves classi- 
fied as of the first kind (Cerveny and Ravindra (10), 

pp. 123-127), in which the head wave is generated in a 
manner shown in Fig. 1.1. The velocity distribution of 
this interface is taken to be Ao > By > A, > By- This as 
the most general distribution resulting in all eight types 
of head waves being formed. The first diagram in Fig. 1.1] 
is for a P wave whose wave front is labelled 1 incident 
upon the interface. Wave fronts which propagate as P waves 
are shown as solid lines while those propagating as S waves 
are depicted as dashed lines. The reflected PP and PS 
(Tinwte) w asmwedll casstransmittedsPPeand PSe(13,2 14) are 
Shown. The four possible head waves which can occur are 
labelled using the numbering system described in Appendix 
A.2. The second diagram shows a similar schematic for an 

S wave (2) incident upon the interface. 


The head wave is a conical wave and in three dimensions 
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Figure 1.1. Schematic picture of incident,,. 


reflected, transmitted and head waves. 
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can be depicted as the frustum of a cone. Therefore, only 
one head wave segment can be produced per ray since the 
spherical wave has been transformed to a conical wave 
front with an infinite radius of curvature in the plane 

of incidence. 

Velocity distributions at the refracting interface 
are vital to head wave formation. Table 1.1 shows the 
possibility of forming any of the eight basic types of 
head waves as a function of velocity distribution. a and 
B represent compressional and shear wave velocities. The 
subscripts refer to the upper (1) or lower (2) medium. 

There is a region in which the impulsive head wave 
always interferes with its corresponding reflected wave. 
This is commonly referred to as the critical region or 
interference zone. It lies immediately beyond the critical 
point r = r*. At the critical point the travel times of 
the head wave and its associated reflected wave are identi- 
cal. The difference in travel times of these waves, 


are lie A 
c - Enands beyond r = r is a positive monotonic- 


reflected 
ally-increasing function of r. If the duration At = I/f 
(f = predominant frequency) of the source pulse is larger 
than this time difference the two waves interfere. If it 
is smaller, both waves exist independently. Hence, the 


end of the interference zone is determined by the expression 
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Table 1.1. Velocity distributions and head wave 
existence. This table displays the possibility of 
formation of any of the eight basic types of head 
waves resulting from an R-wave interacting with any 


of the six types of interfaces. 
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In general it is not possible to find an explicit formula 
for the length of the interference zone and squaeton lz 
must be solved numerically. Appendix A.4 gives more 
details on this. 

A schematic picture of the interference phenomenon 
is shown in Fig 1.2 (from Cerveny and Ravindra (10), p. 168). 
In region I only the reflected PP wave exists. In region 
II the reflected PP and head PPP interfere with each other. 
In region III both waves exist independently of each other. 
Amplitudes of head waves calculated using asymptotic ray 
theory cannot be used in region II due to proximity to the 
singularity r = r*. It is with the total interference wave 
amplitude in this region that this study concerns itself. 

The possibility exists that the receiver may lie 
within the interference zones of the reflected wave and P 
type and S type head wave. This would occur only if the 
distance between first and second critical points was less 
than the length of the P-type head wave's interference 
zone. In this case, all three waves would interfere. 
However, this is a rather artificial situation since it 
requires a large velocity contrast at the refracting inter- 
face (or small difference in P and S velocities in the 
refractor) as well as a very low source frequency. This 
phenomena will not be considered in this study, however 
the formulae for it could be easily obtained from the 
mathematics contained herein. 


It was necessary to calculate the expressions for 
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Figure 1.2. The interference zone. Interference 
zone is denoted by II. The head wave exists only 


beyond the critical point. 
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wave amplitude for all eight head waves listed in Table 1.1. 
To maintain brevity in this publication, only two cases, 
PPP and PSP, will be analyzed in detail as the other six 
cases follow the same formulation. The PPP case has only 
been solved for a multi-layered medium with exactly two P 
ray segments per layer and with only one critical point 
(Cerveny (9)). It is the object of the author to present 
theoretical formulae and numerical evaluations for the 
general problem consisting of any of the eight types of 
head waves with arbitrary velocity distributions at any 
interface. Multiple head waves as described in Chapter 2 
will also be considered. The formulae for interference 
wave amplitudes of all eight cases are presented in 
Appendix A.5. These represent the main new contributions 


by the author. 
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CHAPTER 2 


A WAVE METHOD SOLUTION FOR RAY AMPLITUDE 


eral Introduction to the Exact Wave Method 

The medium under consideration, must consist of m 
elastic, isotropic, homogeneous layers over an elastic 
half-space. The interfaces separating these layers must 
be plane and parallel. At the free surface let there be 
a source of harmonic waves. Also let there be a receiver 
at a given distance, r, from the source and also at the 
free surface. The coordinates in use shall be cylindrical, 
(yr, z, vw) with the origin on the first interface vertic- 
ally below the source. This implies that the source is 
is thesthirckness of the first jayer. 


nike, I ae 0) where h 


18 1 
It will be assumed that the source is symmetrical so that 
the solution is y-independent. The velocities of propaga- 
tion of compressional and shear waves, and density in the 


.th 


S layer will be denoted by Oa 8B ele) GPS LAS eda hs 


June 

respectively. The velocity distribution at the bottom 
interface will be vital in determining what type of head 
waves will arise from the bottoming interface. By consult- 
ing Table 1.1 we find that. the velocity distribution a 
Bat Oe: BF will produce all eight types of head waves. 
In =the ensuing discussion, this 1s the velocity distribu- 


tion which will be taken. 
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The kinematic properties (i.e., travel-time curves) 
of any ray can be totally expressed by the P and S wave 
velocities and the number of P-type and S-type ray segments 
in each layer. The latter quantities are denoted by No j 
and No AC j=1,2,...,m. Figure 2.1-summarizes all the quan- 
tities discussed above. 

It is desirable to derive expressions for horizontal, 
u, and vertical, w, components of displacement due to the 
arrival of an arbitrary ray at the receiver. If the poten- 
tials of compressional waves and shear waves are repre- 


sented by © and ¥ respectively, then the following equations 


are valid in all m layers and the half-space (Cerveny (9)) 
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where © and ¥ satisfy the following wave equations in all 


layers and the half-space 


IMaDur gacucawerckat edie okut 
Greer haf 9 2° af 3 te 
jie 7 ¥ Aan Sa ees 
1 oo (rp 24) + 2S : 2.2 
el a Zz GB, ot 


4 qh xd batonso “ spb theniip ratte! oat 
ar ied sostudsitas 1.4 a a 
svode 

-teamastrod +07 seteeerans evtreb oF whdatath er tL 
ett oF sub asmasaTqetb 46 Btnsnognod yw .footdasy bre 
cuivon sd 42 tovisosh wid Se Gon qodtaldvs ba to TaviaNe, | 
“9494945 ote eovew veae bas gevew (snotezesqmos to zistt 
enoitsups ynfworle} ads nbdt pufevtiosqaey # has # yd begnsa 
({@) yndvxe2) sonqe-tTsd ead bs evever o ffs nt bifsy soa 


:. 
eee le 
i.8 | Ce ) 
| / ove 


ifs nt enotteups svewpmtwollot siz ytatiee ¥ bas @ exenw 
e2sq2-Tlat ot bas etayel 


= ** - 


ees . 


epee ew « 


ee 


18 


source 


receiver 


ert Boel Pu+1 
halfspace ©O 


Bigurcece |. es Nhe: Wed VUMRCeS Chip tiOns. 


Re hid 
rh si pala 


epee ane 


ah =e oer nin CRITE oe iy 
Py pA rs er i es ae a : oa 7 ". 


Le a | 
: Sieh nor 


19 


In this chapter we will only consider those rays 
which have a bottoming reflection of type PP. The PPP and 
PSP head waves will be discussed in reference to their inter- 
ference with the reflected PP wave. The method of calcula- 
tion basically follows Cerveny (6,8,9). The technique has 
been adjusted and considerably expanded to allow an arbitr- 


ary ray. 


2.2 The PP Reflected Wave 

The criterion for this wave to exist is that the 
bottoming reflection be of type PP, i.e., an incident P 
wave and reflected P wave from the interface between layer 
m and the half-space. It is assumed that the source radi- 
ates compressional (P) waves and that their potential 2 


is given by 
rab e | e e 
aye (aak Rg) exp(i k Ro LT we) Zed 


where Ro = (r? + (hy = 2)*)'/2 


w = 2a f = 2 t(predominant source frequency) 
k = w/o, 


By employing the Sommerfeld integral (see Appendix 


A.1) the incident wave potential can be written as. 
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me 1 
Py f Hy (k rv q1) exp(i k[z-h, | (1-95) 1/4) 


=-0O : 


q,(1-q6)7'/? dq, 2.4 


where Ho is a Hankel function of first kind and zeroth order. 
For the duration of these theoretical developments, 
a basic reflection will be defined as the bottoming reflec- 
tion from the interface dividing the mth layer from the 
half-space. 
According to Brekhovskikh (4), the following integral 
expression describes the potential of the wave with a basic 


PP reflection as measured at the surface: 


Raee t <2 1 . — 
os : Ry,(q) o(q) Ho(k r q,) exp(i k B(q))) 


qy(1-a5)"'/* day 2.5 


where R,,(q) is the reflection coefficient of plane waves 


at the mth 


interface. o(q) is the product of all reflection 
and transmission coefficients suffered by the ray in its 


path from source to receiver, except Raq q is defined by 
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The reflection coefficient can be expressed as 


Ry1(a) = Vola) + Vy (a) S(q) + Vo(q) R(q) 


+ V3(q) S(q) R(q) 
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For a full derivation of this consult Appendix A.2. 
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Because of the square roots in this coefficient, the 


Imtegral in 2d BUS ea multi-valued function. 


surface is defined by: 
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arg(v* - ge) le? 


= 1/2 for Qa8 rvs HR oc 65 0 


iF i = 
Ofeaneavbitrary v (say v= IF Oe/Os 72 a./8.47)- 
Employing equations 2.1, we can determine the compon- 


ents of displacement of the wave 


k peel wo t © ] = 
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Bic =P diye 
w= tp t R1,(q9) o(q) Ho (k r qy) 


exp(i k B(q,)) q, 49) 


where Hy denotes a Hankel function of first order and first 
kind. If we assume that we can make a change of integration 
contours for which |q,| > q)> q, © real numbers, for all 
points along the new contour, D, then |k r q,| ker. qy- 
rigker q) Sool ee Chenm ike: q1 | >> 1 and we can employ the 
following asymptotic expansions for the Hankel function 


(see Abramowitz and Stegun (1), p. 364) 
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Hy (x) = (21/2 eyes Gt oe 31) 


with the source and receiver at the same distance from the 


first interface we have 


u = exp(-iw t + U™)(p Ea)! f yy (a) oq) expli k B(aq)) 


(ea ee aH dq, 
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The choice of the contour of integration, D, requires 
some discussion. The original contour (=e.%@) along the 
real axis is unsuitable for three reasons: 


Gh) the integrand oscillates rapidly near the singular 
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DOMRLeS; 

(2) the value of the integrals along this path yields 
displacements for the reflected wave in the broad 
sense (i.e., the sum of the contributions from simple 
reflected wave, head wave, Stonely wave and inhomogen- 
eous head wave). In this investigation interest is 
focused on the contributions of the individual wave 
terms; 

(3) the asymptotic expansions of the Hankel functions 
are invalid along the previous contour. 

The integrand contains four square roots giving it 
its multiplicity; these are stated below with their corres- 


ponding branch points 
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In this investigation only homogeneous head waves of the 

first kind will be considered which eliminates the need 

to study thegeffiect of the last’ two radicals, Peand Q: 
The method of steepest descent is used to calculate 


the integrals. The saddle point is evaluated by solving 


THiSepO Nt, Q] = Xy: is given by numerically solving 


m Nit ereh 45x m Navarh. 2X 
r= J jit ) tT Ls 
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By inspection of 2.8 we see that x, must correspond to some 


angle-related quantity. It is easy to see that Xy is defin- 


able as 
Caras in. 00 
sr dll 2.9 
] Vy 
where oh pe Mike angle of incidence of first ray segment upon 
first interface 
Vi a velocity of propagation of first ray segment. 


This parameter can be transformed to an analogous quantity 
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in any layer by Snell's Law: 


are called the critical angles for the formation of PPP and 
PSP head waves from the bottoming interface. The correspond- 
ing critical parameters in the first layer, xT» 1G can be 


written 


th 


‘ re ake oe wae 
and in the m Taye xX = =300 [0472 xX = a /B The criti 


m m+1° 


cal distance for this ray is the minimal distance from the 
source at which the head wave can appear. Using the quanti- 


ties defined above, we can express these distances as 
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The new integration path consists of several parts 
as shown in Figure 2.2. The portion Dy is a semi-circle 
in the upper half-plane of the complex variable qy° The 


path D is given by the parametric equation 


(lis cE Ue = {Qe ei) + p + exp(-i 1/4); 
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in the real parameter, p. This path passes through the 
saddle point q7 = %] for p = 0. Depending upon the value 
of Xy> and hence upon the epicentral distance, it may be 
necesary to introduce additional contours D, and Do» which 
circumvent the cuts of the radicals S and R. If these cuts 
had not been circumvented the change of variable of inte- 
gration would have resulted in ending up on.a different 
Riemann sheet than we had begun on, making the deformation 
of the path invalid. The branch cuts of the radicals S and 


R are given by the following parametric equations: 
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for p in the real range (0, ~). 

Upon increasing the radius of the semi-circle to 
infinity, the contribution due to it vanishes. The original 
integral is then equal to the integral along D, or along 
D and Dis or along D, D, and D., depending upon the mutual 
position, of x,, Xs xe The following three cases can 


occur: 
(1) GS x < Gy 2] 


In this region only the reflected wave exists, the displace- 


ments of which are denoted by u = y?, Ww = w? 


k 


u? = exp(-i w t + 42)(,* 1/2 { Ry 1(a) ofa) exp(4 k B(q)) 
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(1-q) qy 


(2) xq < xX, < x7 al 


In this region, the epicentral distance lies beyond the 
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branch cut of the radical S. The displacement components 


of the PPP head wave will be denoted by u*, w* and the total 


displacement field by u = u? + u*, w= we + w* where 


Uemome xD cae ett TVs < =e : R1,(q) o(q) exp(i k B(q,)) 


a 
(1-92)71/2 g3/? aq, 


wees exp(-i wt + G(r ea)! f Ryy(a) ofa) exp k B(q,)) 
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The branch of the square roots in the integrals for uP? and 
w? for p = 0 is given by arg(S) = 1/2 and arg(R) = 0 = 
weenie, Ty. 

In the integrals u*~ and w* the arguments of the 
square roots are given by the following for p = O (i.e., 


for q, = x7) 
Mendes U5 UP = 16 


Sulacwaeciffenent sign son each Side ot ithe branch cut. ssee 


the section on the PPP head wave for details. 


* ok 
(3) Xy < Xj <Xy Sq! 


7 


ae 


- or, \ 


4 4 a 


= svertw Sy 4 Ba ake tue? v 


ey Oe 
a. 


(Cypde a» Hake (pds wie MG 


an SNE, SVP =pS pry 


(ippte 2 rh dwadlpia Cs). 4 \" ides Sen! cee 
= 7 - : 

ats a: te : 7 
# Lig > 7 . 


the oy "Oo? ela tepen ada nf efgo* sHuO ad? to tonstd t 
a ; ‘a 
- 0 = (Hote itne Six = (Zhen xd aevtg et 0 = g sot ae 


3] 


In this region the epicentral distance lies beyond the 
second critical point and the contour Do intersects the 
cuts of both radicals S and R. The total displacements 


A * raat - 
are given by u = uo tus tau), wee w? + w + w” where 


u* = exp(-i wo t + Q™)(gap)'/ [ Ryy(a) ofa) expt k B(qq)) 
(leds) mac qe dq, 
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The branch of the square roots in the integrals u?, Wet Or 
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conditions hold for u*, w~ as in the previous case. In the 
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integrals u - w* the arguments of the radicals are given 


by the following: 


The radical R has a different sign on either side of the 
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In this section, displacement components of the PP 
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For this purpose a saddle point method of evaluation will 
be employed. To do this the function B in the exponent of 
the integrand will be expanded in a Taylor series about the 
saddle point hy eS for small values of P. The algebra 


for this is found in Appendix A.3. The expansion is 


B(q,) = B(x,) + i p*(r - r)/2 x3 2.16 
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By using Snell's Law and carefully considering the basic 


geometry of the ray and the medium, it can be shown that 
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which is the travel time curve of the simple reflected wave. 
We introduce the new variable, p, into the integrals in 


2.13 to 2.15 by the parametric equation 2.11. This leads to 
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0 t = 
ur = exp(-f w(t-ty,)) (gag) '/? Jf Ry4(4) ofa) 
exp(-k P=t) p2) gl/2 ap 
2 Xy 
wl = exp(-i w(t-ty3)) (gaz) /? f Ryy(a) ofa) 2.17 
y lye 
exp(-k {t=F) p?) uy p 


where the variables q and q] are functions of the variable 
DaMimaetasinonmaccOraingetoec.th. @lhe Gegionsnear: p = .0 
will have the most significance due to the exponential term 
in the integrand. The functions in the integrand which do 
not vary rapidly with respect to p may be substituted by 
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For brevity, let I denote the integrals in 2.18 


I= f 


Expansion 2.6 is used to express the reflection co- 
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efficient in a more convenient form. Since the Vets are 


smoothly varying functions of p, we can approximate them 
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An approximation can be made to the three integrals in 
equation 2.20 preceded by Vy» V5 and V3. We shall examine 
the approximation in detail for the integral preceded by 
Vy Since the other two follow the same format. Using the 


above substitution for S(p) one arrives at an integral like 
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Then the integral becomes 
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cussed shortly. With € large, one can approximate the 


first radical as 
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for purposes of integration. 

For & to be large compared to unity, k(r-r) is required 
to be large. Since k = 2 w F/o15 it is obvious that the 
condition is valid for high frequencies. From Appendix A.3 


(equations A.3.3 and A.3.4) we see that (r-r) is 
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Since x, < 1, the term (x, - x) will always be positive, 
so that (r-r) is a positive quantity. For the one layer, 

P rays only case, it is clear that r - r = r so that we 
limit ourselves to epicentral distance far from the source; 
far in terms of the wavelength of the source. For the case 
where rays with only P segment are studied, equation 2.16 
shows that r is negative, thus r - r > yr. This is for com- 
pressional velocity increasing with depth. For the general 
situation of mixed P and S ray segments in a multi-layered 
medium, the expressions become too cumbersome to examine 
analytically. The criterion of k(r-r)/2 Xe = m(r=t)/x> r 
>> 1 is tested computationally. Finally, after employing 
this approximation for S(p) and a similar one for R(p), we 


arrive at the following expressions for the radicals: 
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Two cases must now be considered separately. Either 


the epicentral distance is near the first critical point 


~*& 
(x, = x7) or near the second critical point (x, = X1). As 


was discussed in Chapter 1, it is necessary to assume that 


the critical points are sufficiently far apart. 
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"sufficiently" is qualified by the condition that the two 


interference zones do not overlap. 


2.3 Receiver Near the First Critical Point 
Consider bs near the first critical point. In general 
|b. | should be quite large in this region. Expressing 


this mathematically 


This conditicn can only be met for large k. Since k = 

Zar F/a, thasmimpliess that. this result ise Valid only tor 
high frequencies. The termr - r also provides the require- 
ment that the epicentral distances be large. If these 


conditions on b, are satisfied we can set 
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The reflection coefficient Ray can be decomposed as follows: 


Ray = Vo + Vy S + Vo R + V3 Res 
(Vo “ V5 R) - (-V, - V, R) S 
= C - D S 
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With this, equation 2.24 can be written as below if we 
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- D(x) (x = x°)/% (ui (by) ~ 171 goes 


2.4 Receiver Near the Second Critical Point 


In this region we assume 
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As in the previous case this condition requires a high 


frequency for fulfillment. With these conditions we can set 
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Applying these approximations to equation 2.22 yields 
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The reflection coefficient can also be decomposed as 


Ray = Vo + Vy SH Vo Reet V3 Sue OS) 
= (VG - Vy S) - CNG - V3 Spar 
= Cl - Dl R 


Ry, (x) = C1(x) - ASG @ INO, 


We again employ Uy as defined above, with the result 
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At this point we shall summarize by stating the 
horizontal, u?, and vertical we, components of displacements 
for the simple reflected wave. We substitute 2.25 and 2.26 
into 2.18. L is the geometric spreading factor. In Appendix 
A.3 it is demonstrated that it is very nearly equivalent 
to that derived by asymptotic ray theory. In this develop- 


ment we define L as 
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Nearathe second critical point: 
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2.5 The PPP Head Wave 

In this section the expression for displacements of 
a head wave, propagating along the interface between the 
meh layer and the half-space (layer m+ 1), as recorded 
at the surface, will be discussed. The expressions to be 
evaluated are given by 2.14. We introduce a new integra- 
tion parameter, p, defined by 2.12(a), into 2.14. The 
values of p, on both sides of the branch cut will be in the 
range (0, ~). The value of S(x) will have different signs 


on either side of the path. (S(X)) yex(p) can be extended 


using 2.12(a) 
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where 
0 on left bank 


TT on right bank 
i * 
arg(p + 2 e! es (1 - eee 7/0 fOveep. =a0 


In the same fashion R(x) is expanded 
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B(q,) is expanded about xy The derivation is in 


Appendix A.3 
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By considering the geometry of the medium and the ray 
travelling through it, the formula for the travel-time 


curve of the PPP head wave can be set down as 
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and 


We now introduce the concept of head wave multipli- 
city. Figure 2.3a depicts a ray which has M PP-type 
reflections from the bottoming interface. Any one (and 
only one) of these reflections can form a head wave which 
will be returned to the surface. The reason that only one 
reflection can form a head wave follows from wave front 
considerations. Head waves arise from a wave front of 
finite radius in the plane of incidence impinging upon a 
plane interface. The head wave is a conical wave in three 
dimensions. Hence it has an infinite radius of curvature 
in the plane of incidence. For this reason it cannot form 
another head wave. The ray depicted in Figure 2.3b is a 
physica bhysimpossible=sT tuations = This=is often mistakenly 
conceived in the literature (see Hales and Nation (13), 

p. 532). To properly evaluate the head wave displacement 
we must consider (Rien rather than Ray in the integrals 
Cola. 

Consider any two rays which travel the same ray path, 
but with one difference. The two rays form head waves 
from different bottoming reflections. These two rays have 
identical travel time curves (hence the rays are termed 
KINEMATIC ANALOGS). They suffer the same reflections and 


transmissions and their amplitude-distance curves are identi- 
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Figure 2.3. Head wave multiplicity 
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cal. In this case they are further termed DYNAMIC ANALOGS. 
For more details on kinematic and dynamic analogs refer to 
Hron (17). Within the critical region for the reflected 
wave-head wave pair each member of the group of dynamic 
analogs interferes equally with the reflected wave. 
Physical reasoning leads one to the conclusion that the 
total displacement of this group should be M times the dis- 
placement of an individual member. We now give a more 
rigorous mathematical derivation of this. 
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Using equation 2.6, the M power of Ry] can be 


written as 


[R, (a) ]" = [(Vgla) + Vola) R(q)) 


(-V,(q) - V3(q) R(q))S(q) 1" 
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Near the first critical point the radical S is very nearly 
zero. Thus the second term in brackets is very much less 
than one. The expression can then be expanded in a bi- 


nomial series omitting the third and higher terms: 
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The functions V. do notecontain thesradicates or Rw and 
hence are slowly varying functions. The saddle point 
approximation allows one to replace these values in the 
integrand by their value at p = 0. The same approximation 
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Some further approximations can be made on equation 
Co Oem | NCCn Om Neawetnhemfinste chi ticals points mite Us 


justified to approximate the reflection coefficient as 


C does not, however, contain the radical S and therefore 


has a finite derivative near the critical point. cM-1 can 
then be replaced by its value at p = 0, namely cM-1 = 
M-1 


Ri] . This term represents all the bottoming reflections 
of type PP which do not form a head wave. For the duration 
of this discussion the term Reg i(x) will be included in 


the product AG). Thus, equation 2.35 can be rewritten as 
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Next expand the individual terms and integrate, remembering 
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that the arguments of the term p in the expansion of 


S(q) has a different sign on either bank of the cut 
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In the previous expressions we have used q as a vari- 
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able, but it was understood that q is a function of p in 
a manner according to equation 2.12a. Since the V's do 
not contain radicals R and S, they can be substituted by 
their value at p = 0. Using the full expansions of R and 


S as given by equations 2.29 and 2.30 we have 
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We can rewrite the displacements as follows: 
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in the exact formula 2.38. 


2.6 The PSP Head Wave 

The PSP head wave exists only for epicentral distances 
greater than the second critical point or analogously, for 
ay ie To calculate the displacements for this wave, 
we must evaluate the integrals in equation 2.15 where dD. 
is the contour which circumvents the branch cut given para- 
metrically by 
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The method of evaluation is very similar to that used in 
the previous section so that only the major highlights and 
final results will be discussed here. 
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where for p = 0, both radicals have argument of zero. 
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Upon substituting these expressions into equations 2.15, 
and by allowing slowly varying functions of qy to assume 
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The multiplicity factor, M, has the same significance 
here as in the case for PPP head waves. The term (Ree q) oe 
is binomially expanded and apprxomations similar to those 
for PPP head waves are made. Expansion 2.6 is employed 
for the reflection coefficient and the V's are assigned 
GheLoey a luemat ps -2 0s qi = oe o,/8.47)+ Rs which has the 
Same modulus on either bank of the cut, has a different 
argument by mt on either bank. For this reason Vy drops as 
a Peedi of integration and V5 and V3 survive integration. 
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the integral results from the region where p is near zero. 
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stitutions for the second critical point made. If the two 


critical points are adequately separated, or in other words 
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These equations again have a singularity at the critical 


pOInNt Whicheis not presentein’ 2,45. 


2.7 The Interference Wave 

In Chapter 1, the relationship between head waves and 
simple reflected waves within the interference zone was 
discussed. Thus, we must consider the summed total response 
due to the head wave interfering with the reflected wave. 
The previous discussions have determined that the total 
response of multiply-reflected head wave is M times the 
response due to a single head wave. With this in mind, the 
displacements of the interference waves can be written as 
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Only the calculation for the PP-PPP case will be 
given in detail since the PP-PSP case follows the same 
procedure. The displacements referred to in the equations 


aboveware defined in equations ¢.2/, 2.28, 2.38 and 2.45. 
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Since the receiver lies in the vicinity of the first 


critical point we can make the following approximations: 
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The variable in G in this case has been left as by to ensure 
more accuracy. 
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Equation 2.50 is valid at all epicentral distances. 
We will note in the following chapter that the complex 
valued function G(x) tends to zero for large values of x. 
Hence, the second term decreases in significance with 
respect to the reflection coefficient Ray for epicentral 
distances far from the critical point. Since the second 
term is inversely dependent on source frequency, f, it is 
obvious that for the high frequency case the corrective 


term nearly vanishes leaving 
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using asymptotic ray theory for a high frequency. For this 


case the head wave amplitude is negligible. The results 
of computations presented in the next chapter bear this 


out quite well. 
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wave in the neighbourhood of the critical point can be 


described by 
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where aM and abr are the modulii of the complex amplitudes 
in the horizontal and vertical direction. Z represents 
the argument of the complex number. To facilitate the use 
of one common term it proves expedient to introduce the 
concept of total surface incident amplitude, expressible 
as a complex number. For purposes of this discussion refer 
tosrFigure 2.4. 

The total wave system generated by a wave incident 
upon the free surface is illustrated in Figure 2.4(a). 
What the receiver (R) records is actually the interference 
of the incident wave and the reflected P and S waves. To 
calculate the effect of the earth's surface on the recorded 
wave, the incident amplitude must be multiplied by an 
appropriate surface conversion coefficient. These are 
listed in detail in Hron (15). The surface conversion 
coefficients also resolve the recorded total wave into 
horizontal and vertical components. 

For the remainder of this thesis only the total inci- 
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sidered. For the sake of brevity, we have only discussed 
two of the eight head waves listed in Table 1.1. A similar 
evaluation was performed for the remainder. Upon complet- 
ing these non-trivial computations, the complex amplitudes 
of all eight interference waves can be described by the 
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yp for r< r 
Jes 
( Femi aWe ECtiaves t 1/2 ton p> 
apn kas = 


tag oting = travel times of simply reflected wave (or 
reflected wave in the narrow sense of the 


word) and head wave 


Dike = head wave coefficient. They are derivable 
from the reflection coefficients R; Q 
(see Cerveny and Ravindra (10), p. 131f 


or Hron (15)) 
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which is incident upon bottoming inter- 


face and which forms the head wave. 


The parameter x as defined above is the one which was 
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used in all numerical evaluations. This use of a parameter 
will permit the same branch cuts in the complex plane for 
incident S waves as incident P waves. 

The next chapter provides numerical evaluation of the 
quantities calculated here. In most instances a close 
resemblance between this theory and asymptotic ray theory 


outside the critical region can be noted. 
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CHAPTER 3 


NUMERICAL RESULTS 


3.1 The Weber Function 

The function G(y) as defined by equation 2.49 can 
be expressed in terms of an appropriate Weber function (or 
alternatively named a Parabolic Cylinder Function). The 
following integral definition is essential (see Magnus 


and Oberhettinger (19), p. 92) 
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ee a D(z) 
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Using this definition, G(y) may be rewritten 


‘ see 
exp(Q® - EY) dy o(y(i-1)) - 4 29/4 yl/? 3, 


where G(y) is a complex valued function of the real para- 
meter y, -~> < y < ™, 
Kireyeva and Karpov (18) have tabulated the proper- 


ties and numerical values of the Weber functions. These 
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functions arise when one solves the Helmholtz equation 


in parabolic coordinates by a separation of variables techni- 


que. This results in Weber's equation 


2 


~ =) D.(z) = 0 
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Kireyeva and Karpov have solved this equation for z = 
y(1+i) and real t. The real and imaginary parts of 
D.(y(1+i)) = u,(y) + j v,(y) solve the following system of 


equations: 
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ot/2 (2 n)'/2 


v,(0) = 0 v;,(0) = « F(-t72)} 


By using two important properties of the Weber function 


(Kireyeva and Karpov (18), p. x). 


D,(z) De (z*) 


D.(z) = (-1)* D,(-z) 

we can find the solution for Dy poly(i-1)). To solve this 
problem numerically, one may utilize SUBROUTINE DASCRU 

from the International Mathematical and Statistical Library 
(IMSL) or any similar subroutine in standard computer 
libraries. This routine provides a numerical solution to 

a given system of first order ordinary differential equa- 


tions. By introducing 


dmu (y) dv (y) 
ue a is 
or as = ay y) —————— = 


into equation 3.2, the system is transformed to first order. 
By numerical tests two important results were 
obtained: the step size in the variable of the Weber func- 
tion and the range of the variable. By checking the com- 
puted results with the results of Kireyeva and Karpov, it 


was determined that a step size of .005 for y would facili- 
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tate linear interpolation in the functional values. For 
most seismological situations, the parameter y, defined 
in equation 2.48, within the range (-5,10) will sufficiently 
cover the interference zone and epicentral distances near 
itee Of fcourse, y = 0 at the@cwmitical point. 

A FORTRAN program has been created to compute the 
Weber function and write the value (unformatted) onto a 
disk file for future retrieval and interpolation in the 
construction of amplitude-distance curves. Using the un- 
formatted write facility of the Michigan Terminal System 
(MTS), up to 248 bytes per line of disk file can be utilized. 
The program writes 30 pairs of single precision real numbers 
(i.e., representing real and imaginary parts of 30 complex 
valued numbers) onto one line. The last complex number on 
line j is repeated as the first complex number on line 
j + 1. This is to facilitate interpolation between any two 
values of a Weber function by reading only one line of the 
disk file. Another FORTRAN program was created to read 
the Weber function file and determine the functional value 
for any parameter. Given a certain value of the parameter, 
y, the program determines and reads the line of the file 
which contains functional values for parameters closest 
to y. Subsequently it selects and interpolates between the 
appropriate functional values. The derived Weber func- 
tional value at y is then adjusted by the proper quantities 
to form G(y) as defined in equation 3.1. 


Figure 3.1 is a computer generated plot of G(y) pro- 
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duced on the CalComp plotter. Before the critical point 

(y < 0) the real and imaginary parts of G(y) increase-in 

a continuously monotonic fashion. Beyond the critical 
point (y > 0) they decrease rapidly and oscillate about 
zero. For y < -3 the real and imaginary parts of G(y) 
become very insignificant. Consequently, the second term 
invequation’2°95lwill®become trivial ¥compared “to "the reflec- 
tion coefficient. The equation degenerates to the value 
for the reflected wave derived using asymptotic ray theory. 
Alternatively, for y > 7 the function attains very low 
values and again equation 2.51 tends toward the asymptotic 
ray theory formula. It is also advantageous to take cogni- 
zance of the fact that the function attains maximum 
modulus at y = 0 (characterizing a critical point) and 

thdt Tt “oscillates beyond Mt#e Thus;y one can “infer ‘that 

the amplitude-distance curves will behave in a similar 


fashion. This will be displayed in the following sections. 


3.2 Computed Dynamic Properties 

A FORTRAN program has been written to calculate 
modulus and argument of the amplitude of interference 
waves. The results can be subsequently plotted by a Cal 
Comp plotter. The basic formula to be evaluated for many 
epicentral distances is equation 2.51. 

Figures 3.2-3.4 are a series of diagrams displaying 
the logarithm of the modulus (base 10) and argument versus 


distance from the harmonic source. In all of these figures 
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the source frequency (FO) is given at the top as well as 
the type of interference wave. The interference modulus 
curve is marked in a heavy solid line, whereas the head 
wave (which can only appear beyond the critical point R*) 
has its modulus shown in the light line. This is the head 
wave modulus as calculated by asymptotic ray theory. The 
simply reflected wave modulus as computed by the ray theory 
is shown as a dashed line. The end of the interference 
zone is indicated by RI. 

The argument plot, lying below the modulus plot 
gives the argument (expressed as a value between 0 and 2 7m) 
for the interference wave (X), the simply reflected wave 
(Ojmandmthesheadewavemcassolidsdotwat ther criticalspoint). 
Since the ray parameter for head waves and head wave co- 
efficient are independent of epicentral distance, the argu- 
ment is the same for all distances. The elastic and spa- 
tial parameters of the medium are given at the bottom of 
the diagram. Velocities are expressed in km./sec., densi- 
ties as specific gravity and layer thickness in km. 

For the remainder of this discussion the term DYNAMIC 
PROPERTIES will refer to the modulus and the argument of 


a particular type of wave. 


3.3 Dynamic Properties and Frequency 
Figure 3.2 shows dynamic properties for the interfer- 
ence of a PP reflected wave and a PPP head wave for four 


frequencies. The asymptotic curve for the simply reflected 
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Figure 3.2. Dynamic properties and frequency. A series 
of plots of dynamic properties is displayed for an 
interference PP-PPP wave system for frequencies of 
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ray has a sharp peak at the critical point and continuously 
decreases beyond this. The interference curve is smooth 
in the vicinity of the critical point and attains a maximum 
somewhat beyond the critical point. By close examination 
of Figure 3.2, it can be seen that the position of the maxi- ~ 
mum shifts closer to the critical point as the source fre- 
quency increases. This result was also predicted by 
Cerveny (4) (pp. 215-219). At greater distances, the inter- 
ference curve decreases as it oscillates about the ray 
theoretical curve. With higher frequency the width and 
height of the oscillations in the interference wave curve 
decrease. The argument of the interference wave also 
oscillates about that of the simply-reflected wave. How- 
ever, for the higher frequencies these oscillations become 
nearly invisible. 

The length of the interference zone is inversely 
proportional to frequency. The formulae describing the 
length of the interference zone are discussed in Appendix 


Aisi. 


3.4 Dynamic Properties and Multi-Layered Media 

The application of this theory to a multi-layered 
medium is demonstrated in this section. Let us consider 
a model consisting of two layers over a half-space. The 
model chosen consists of the first two layers of the 
Alberta model described in Hron and Kanasewich (16). The 


velocity distribution at the bottoming interface is such 
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as to allow the formation of all eight types of head waves. 
In all examples shown the simply reflected ray consisted 

of four segments. Two of these were always P-type in the 
first layer. The other two were in layer two and of the 
type stated on each plot. The source frequency is constant 
dium cao sie GtZe lea llecasess. 

It can be readily seen that the fit between interfer- 
ence amplitude and simple reflected amplitude grows better 
with increasing distance from the critical point. Upon 
examining the PP-PPP case it is clear that the simple 
reflected amplitude curve has little resemblance in general 
Shape to that of Figure 3.2. Two reasons can be presented 
for this. The earth model used in constructing Figure 3.3 
has a second critical point which causes a peak near this 
point. Secondly, when computing ray amplitude the trans- 
mission coefficients at the interface between layers 1 and 
2 tend to smooth out the sharp peak which appears at the 
finsivecmicical point. 

The PP-PSP case also requires some discussion. This 
head wave has the largest critical distance of all eight 
types. At larger epicentral distances the spreading fac- 
tor, L, in equation 2.51 becomes more significant than at 
small ones. Thus, for the relatively large epicentral 
distances considered in this case the spreading factor has 
become large enough to damp out any oscillation present in 
the amplitude distance curve. A numerical experiment was 


performed to calculate the value of the amplitude without 
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Figure 3.3. Dynamic properties and multi-layered media. 
A series of plots of dynamic properties in a two layered 
medium for all eight head waves. Figure 3.3c shows the 
case PP-PSP with the effect of geometrical spreading 

removed to display the oscillatory nature of the modulus 


curve. 
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including the effect of the spreading factor, L. The 


resultssare shown in Figurer3.cc. 


3.5 Dynamic Properties and Head Wave Multiplicity 

Figure 3.4 consists of four plots of dynamic properties 
fOmmaehPPahead waves wi theamuditipricity M =mlecs3. deem For all 
cases the source frequency is 15 Hertz. The lines sketched 
On each graph have a somewhat different definition than 
in the previous sets. The solid line running over all epi- 
central distances up to the end of the interference zone 
is the interference wave amplitude calculated using formu- 
lae derived in this thesis. This line becomes a line of 
long dashes beyond this point. This is to signify a cer- 
tain inaccuracy in the formulae which will be discussed 
later. The sclid line beginning beyond the interference 
zone represents the head wave amplitude as derived from 
asymptotic ray theory. The dashed line is split into two 
segments each having different significance. Before the 
end of the interference zone the dashed line represents the 
simply reflected PP wave amplitude computed from asymptotic 
ray theory. Beyond the end of the interference zone it 
represents the sum of the simply reflected and head waves. 
To produce this sum one considers the reflected wave dis- 
placement exp(-i w(t - t__)) > Ap and the head wave dis- 


Pp 


placement exp(-i w(t - t )) Ay where Ap and A, are complex 


ppp 
numbers representing the dynamic properties of the wave. 


These two displacements are summed as 
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Figure 3.4. Dynamic properties and head wave multipli- 
city. A PP-PPP interference wave has its dynamic 
properties plotted for multiplicity.) M;.equalleiO lsc. 


3, 4. Source frequency is constant at 15 Hertz. 
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The subscripts r and i refer to real and imaginary parts, 


respectively. This can be adjusted to 


=1he03. (tt = tien) 


a) PP [(A COSa mes ca Aus S1n @ At) 


: i(Ap, a? ne COS Ne Wwe Sy eh eh 
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7 PRT PAT MRAte= CHle_ we 
7 LAI pp ppp 
Modulus = |A| Argument = tan7! ATs 


On the argument plot the x's again represent the 
interference wave. The open circles denote argument of 
the reflected wave amplitude up to the end of the inter- 
ference zone and the argument of the total asymptotic ampli- 
tude beyond this point. dhe solid circle at the critical 
point gives the argument of the head wave amplitude. 

Several points should be noted about these diagram: 
(is) With increasing head wave multiplicity the upward 

shift of the interference wave amplitude from the 


reflected wave amplitude increases. This is what 
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one would expect on the basis of physical reasoning 
alone. This increasing shift upwards is due to the 
increasing number of head waves with the same ampli- 
tude interfering with the reflected wave. 

(2) With increasing multiplicity the ray theoretical 
total amplitude curve tends to behave more like the 
head wave curve than the interference curve. This 
feature goes hand in hand with the fact that the head 
wave amplitude increases with increasing multiplicity. 

(3) With increasing multiplicity the phase shift seems 
to have a poorer fit between interference and asymp- 
totic ray values. The reasons for this will be dis- 
cussed in point (4). 

(4) The pronounced dip appearing in the interference 
wave amplitude curve can be explained by the influ- 
ence of the free surface coefficient. This coefficient 
undergoes a phase shift where it passes through zero 
modulus (F. Hron, personal communication). With an 
increasing number of reflections from the free sur- 
face this effect becomes exaggerated. The interfer- 
ence curve has been drawn in long dashes beyond the 
interference zone due to this inaccuracy. A fuller 
discussionrotethis#isenecessanry: 

We recall from Chapter 2 that the generalized dynamic 


property formula was of the form 
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where R, g represented the reflection coefficient, H, Ro 
the correction due to the head wave in the critical region, 
L the geometric spreading and o the product of all inter- 
face coefficients encountered by the ray in its path from 
source to receiver. We also recall that o(q,) was removed 
from the integrand and q4 assigned its value at the saddle 
point, namely Xa: Near the saddle point this approximation 
is valid, but with increasing distance its validity falters. 
Thus when the free surface reflection coefficient passes 
through zero it totally suppresses the term Bs k 2 in the 
above formula. 

However, for the purpose for which this study was 
undertaken (dynamic properties in the critical region) the 
results are only necessary within the interference zone. 
Should the formulae be used for the calculation of synthe- 
tic seismograms, their use should be restricted to the 


critical region. Outside that region the asymptotic ray 


values should be taken. 


3.6 Reflected Pulses Near the Critical Point 

On the basis of formulae presented in this study, 
synthetic seismograms can be produced for interference 
waves near the critical point. At points far from the 
critical point, one may utilize the amplitude formulae 
from asymptotic ray theory. The model chosen for these 
computations has been deliberately simplified for purposes 


of display. Figure 3.5 shows the source pulse for this 
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Figure m.5.° Source pulse for synthetic 


seismograms. 
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calculation. The synthetic seismograms for horizontal 

and vertical components of an interference PP-PPP wave 

are given in Figure 3.6. The earth model is the Manitoba 

model (one layer over a half-space) as given in Hron and 

Kanasewich (16). The critical distance (Ro) and the end 

of the interference zone (Ry) are marked on the plot. For 

continuity, the exact formulae were used from .8 times the 

critical distance to the end of the interference zone plus 

.1 times the length of the interference zone (approximately 

from 55 to 126 km.). Outside this range the ray theore-- 

tical expressions for amplitude were used. 
Several conclusions can be drawn from Figure 3.6: 

Gh) The form of the reflected pulse changes little on 
passing through the critical point. The pulse does 
suffer come distortion due to the reflection coeffi- 
cient becoming a complex number beyond the critical 
point. 

(2) The total separation of the head wave from the reflected 
wave occurs only after the end of the interference 
zone. 

(3) The travel-time curves of the head and reflected 
waves are tangential at the critical point. 

(4) Beyond the interference zone the head wave amplitude 
is very insignificant compared to the simply reflected 
wave. The head wave amplitude could be enhanced 
(relative to the simply reflected wave amplitude) by 


an increase in multiplicity. With an increase in 
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multiplicity the head wave amplitude may becomes 

non-trivial next to the simply-reflected wave. 

(5) The head wave pulse is distorted by an argument of 
1™/2 with respect to the source pulse. Figure 3.5 
displays the source pulse as a sinusoidal pulse in 
time, while Figure 3.6 exhibits the head wave as a 
negative cosine in time. In terms of asymptotic 
theory this is explained by the factor of 1/i in the 
head wave amplitude formula. 

It was shown in previous sections that the head wave 
amplitude is inversely proportional to frequency. With an 
increasing frequency the head wave amplitude will become 
such a low value outside the interference zone that its 
contribution to the seismogram will become insignificant. 

Figure 3.6 indicates that there is a very good tran- 
sition between the exact formulae for dynamic properties 
and the asymptotic ones. This again provides us with a 


check upon our derivations. 
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CHAPTER 4 


CONCLUSION 


This thesis presents a theoretical study of the 
dynamic properties of interference head waves near the 
critical point. The results are valid for arbitrary velo- 
city distributions along the vertical in the bottoming 
ave All rays, including those with multiple bottoming 
reflections are suitable to have these calculations per- 
formed on them. 

Chapter 1 gives an introduction to the study of 
reflected and head wave systems. An extensive, although 
not exhaustive, bibliography is given. 

Chapter 2 states the actual theoretical solution 
of the problem. The integral expressions for the reflected 
wave in the broad sense are evaluated along special com- 
plex-valued integration paths. The contours of integration 
chosen are a modification of the saddle point contour, 
allowing certain quantities to be approximated as con- 
stants. Chapter 2 also demonstrates the similarity 
between this method of calculation and asymptotic ray 
theory. This resemblance becomes stronger for high fre- 
quencies or for distances greatly removed from the critical 
Dou nus 


In Chapter 3 the expressions derived in the previous 
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chapter are cast in a form suitable for numerical computa- 
tion. In particular, attention is focused on the Weber 
function, Di se: A series of numerical results are given 
in which dynamic properties of interference waves are 
plotted against epicentral distance. The quantities varied 
include source frequency, number of layers in the over- 
burden, and head wave multiplicity. Also shown in these 
sections are the dynamic properties as computed by asymp- 
totic ray theory. The similarities between these results 
are discussed. Finally, Chapter 3 provides a set of 
synthetic seismograms displaying interference wave char- 
acteristics. 

This thesis has presented new and original material 

in the following areas: 

(1) Formulae have been derived for all eight possible 
types of interference head waves. These formu- 
lae are valid at all epicentral distances, in 
particular the critical region. The formulae 
tend to asymptotic ray theory values for high 
frequency or great distance from the critical 
point. 

(2) The works of Cerveny have been generalized to 
the cases of multi-layered media, arbitrary 
ray paths and any of the eight possible types 
of interference waves. Related to arbitrary 
ray paths are multiple head waves, the solution 


of which is also presented in this study. 
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(3) Numerical solution and tabulation of the Weber 
function Dis;2° 
There are several areas of this problem that require 
further attention. Some of the possibilities are: 
(1) Exact numerical integration of the integral 
expressions for the reflected wave system could 
be performed. Although this is computationally 
expensive, it would provide an excellent check 
upon both the asymptotic theory and the one 
presented here. 
(2) A similar approach to this one should be 
attempted for an anisotropic medium. The major 
difficulty here is not so much in solution, but 
rather in formulation of the integral expressions. 
(3) The dynamic properties should be studied for a 
group of kinematic analogs rather than for a 
single ray. These results could then be 
employed in the computation of synthetic seismo- 


grams. 
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APPENDIX A.1 


THE SOMMERFELD INTEGRAL 
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The Sommerfeld integral is defined in Ewing et al 
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Let d,(pr) = w2tH6") (pr) + KH?) (pr) J where Ht) and H{2) 


are Hankel functions of the first and second kind. Refer 


to Abramowitz and Stegun (1). Doing this we obtain the 


sum of two integrals: 
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In the second integral replace p by -p and utilize the fact 
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Two integrals are obtained with identical integrands but 


with limits (0, ~) and (-~, 0). Combining these yields 


Substitute p = q4 k and we arrive at 
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This is the expression for the incident wave potential. 
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APPENDIX A.2 


TRANSMISSION AND REFLECTION COEFFICIENTS 
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All reflection and transmission coefficients which 
were used in these derivations and in numerical evaluations 
are obtained from the work of Hron (15). An in-depth 
account of how the V's are derived from the PP reflection 
coefficient will be given. The other three coefficients 
follow a similar routine. 


Hron (15) makes the following preliminary definitions: 
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plications. 
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Head wave coefficients can be derived from these 
coefficients by using a simple relationship involving the 
reflection coefficient. Head wave coefficients yield the 
relationship between the amplitude of an incident wave and 
the resultant amplitude of the head wave arising from its 
critical refraction. In general any of the four basic 


reflection coefficients can be written as 


where H is either S (when considering a P type head wave) 
or R (when considering an S type head wave). The terms 
HG yell whe - M are in general complex and do not contain the 
radical iH. 

Head wave coefficients can be obtained from reflec- 


tion coefficients in the following fashion (Cerveny and 


Ravindra (10)) 
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the particular type of head wave under consideration. 

The numbering system used to identify the various 
reflection and head wave coefficients requires some explana- 
tion. It is necessary to label the four velocities at the 
basic reflecting interface. These are labelled as ails 
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" melo? Baty te To determine the number of the 


reflection or head wave coefficient one considers the velo- 
city at which each ray segment propagates. Thus a PP 
reflection requires the coefficient labelled 11 (Rj). As 
well a PPP (PSP) head wave has a head wave coefficient 


labelled 131 (141). The PPP (PSP) head wave coefficient 
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THE EXPANSION OF B(q,) 
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B(q,) is defined in Section 2.2 as 
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A new contour of integration and a new variable of integra- 


tion, p, are introduced via 
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The function B(q,) can be expanded in a Taylor series at 


the point 


X41 = Sie 04 SS ore 


where Xy is the root of equation 2.8. The series can be 
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Both derivatives in the expansion are to be evaluated at 
the saddle point, ly te hee The first derivative when 
evaluated at this point gives zero results since this is 


the definition of the saddle point. We are then left with 
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where Po is the value of p at the saddle point. Obviously, 
Po = 0. Substituting these back into the expansion one 
obtains 
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One then multiplies the term (x, - x) into the brackets. 


After doing this it proves expedient to add and subtract 
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Equation 2.8 is then employed. This equation gives the 
relationship between xy and epicentral distance r. As 


well, the quantity r is introduced as 
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The expansion of B(q,) along the contour used to 
evaluate head wave terms will have some vital differences. 
Let us only consider the P-type head wave term. The S- 
type head wave term is similar. The contour in this case 


circumvents the cut given by 
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Thus, the first term in A.3.4 should read B(x>) rather 
than B(x i As well, the second term should read 
* 
Wee xy o oafwwas just ¥ with x, substituted 


"ppp ppp 
for Xy° Going back to equation A.3.2 we shall show that 
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the second term does not vanish for this case. Evaluating, 
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According to equation 2.10 these last two summations yield 
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We now demonstrate the close resemblance between the 
geometrical spreading factor derived in this thesis, L, 
with that from asymptotic ray theory, Lp: The formula 
presented early in Chapter 1 showed the geometrical spread- 


ing from ray theory to be 
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Some adjustments should be made to this formula to make it 
consistent with the notation used in this work. First, 

the summation should take place over m layers rather than k 
ray segments. The number of P and S ray segments per layer, 
N. . and N. ; are introduced. Second, the ray parameters, 


pi 
x, are used rather than cos 0.. These changes enable one 
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to express asymptotic ray theoretical spreading as 
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where Vy represents the velocity of propagation of the 
first ray segment. The geometrical spreading factor 


derived here is expressed as 


Employing A.3.3 for (r-r)'/* and equation 2.8 for r we 


arrive at 
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It is known that at the critical point a head wave 
forms in addition to the reflected wave in the narrow 
sense of the word. Up to a certain distance the waves 
interfere. The length of this interference zone or criti- 
cal region has been ably discussed by Cerveny (7,9) and 
Cerveny and Ravindra (10). The formulae presented in this 
appendix are derived in these works and are presented here 
for completeness. | 

The end of the interference zone is defined as the 
distance beyond the critical point at which the difference 
in arrival times of simply reflected and head waves equals 


the duration of the source pulse. This can be stated as 
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For a multi-layered medium each layer is treated separately 
for each type of head wave that can possibly exist at the 
base of that layer. Two types of head waves must be con- 
sidered; symmetric and non-symmetric. A symmetric head 
wave is one in which the incident ray segment is of the 
same type (P or S) as the reflected ray segment. These 

are the head waves PPP, PSP, SPS, SSS. Conversely, the 
non-symmetric head waves are PPS, PSS, SPP, SSP. For 
symmetric head waves, an analytic expression can be found 
for the length of the interference zone. According to 


Cerveny and Ravindra (10) this length can be expressed as 
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V; represents the velocity of propagation of both incident 
and reflected ray segments of the head wave. Vy is the 
head wave velocity along the bottoming interface. The 
frequency and layer thickness are given by f and h respec- 
tively. 

This formula shows that for increasing frequency f 
the length of the interference zone decreases. 

For non-symmetric head waves we consider the travel- 


times of the reflected and head waves through the layer as 
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ence in travel-times may now be solved numerically for the 
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the ray entered the layer. Subtracting this distance from 
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of the interference zone. ri ane is expressible as 
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For all interference zone calculations in this thesis 
the equations for the zone length of non-symmetric waves 
were solved using a computer library subroutine. The actual 
solution was performed by interval halving. The subroutine 


produces accuracy greater than Haar 
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In this Appendix, we give the explicit formulae for 
calculating the total incident amplitude of any of the 
eight interference waves. The numbering system of the 


waves is explained in Appendix A.1. 
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All quantities have been defined in equation 2.51. 
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